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Abstract 

Collisions of polar state molecules at ultralow energies are considered, within a model that 
accounts for long-range dipole-dipole interactions, plus rotation of the molecules. We predict a 
substantial suppression of dipole-driven inelastic collisions at high values of the applied electric 
field, namely, field values of several times Be/n- Here is the rotational constant, and fj, is 
the electric dipole moment of molecules. The sudden large drop in the inelastic cross section is 
attributed to the onset of degeneracy between molecular rotational levels, which dramatically alters 
the scattering Hamiltonian. As a result of the large ratio of elastic to inelastic collision rates, we 
predict that evaporative cooling may be feasible for state molecules in weak- field-seeking states, 
provided a large bias electric field is present. 
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I. INTRODUCTION 



The paramount goal of physics of ultracold temperatures is to control and manipulate the 
quantum world. Polar molecules bring new challenges and hopes to this field(see the review 
in ^). Since 1998 P], several experimental groups have been restraining polar molecules in 
order to get colder and denser samples. The difficulties on this road speak for themselves: 
a whole variety of novel experimental techniques were developed for this purpose But 
in spite of very intensive experimental research, the production of ultracold molecules still 
poses a significant challenge. Cold and dense samples would allow one to control two- 
body J IJand many-body systems of polar particles 0,y,l2,3,y,Hlllll2,|l3,14,|l3, 
although so far none of these goals has been experimentally realized. The main obstacle 
to these achievements is loss of trapped molecules via inelastic collisions, or else Majorana 
transitions [l^. In this paper we study just one aspect of physics of cold polar molecules, 
namely, the possibility to have molecules in a weak-field-seeking state with sufficiently low 
inelastic collision rates to allow evaporative cooling. 

We have previously considered the electrostatic trapping of polar H-state molecules of 



both bosonic (OH) anc 

r 

respect to collisions [17, 



ermionic (OD) symmetry from the point of view of stability with 



18|. As electrostatic trapping requires molecules to be in a weak- 



electric-field-seeking state, collisions involving the strong and anisotropic dipole-dipole inter- 
action between molecules may drive the molecules into unfavorable lower-energy strong-field 
seeking states, leading to unacceptably high trap loss and heating. For bosonic 11- state OH 
molecules we have found that the elastic rate can be much larger than the inelastic rate only 
for quite large field values. As the first excited rotational level of OH lies 847^^ above the 
ground state, inelastic rates are defined mostly by A- doubling and a hyperfine splitting. In 
general one cannot yet exclude the possibility of finding molecules whose hyperfine structure 
permits evaporative cooling, but such a candidate has not yet been identified. 

Polar fermions have a potentially important advantage for electrostatic trapping, namely, 
low inelastic rates at cold temperatures. The state-changing collisions of dipolar fermionic 



molecules were discussed in 
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20|. Based on the well-known Wigner threshold laws for 



dipole- dipole interactions it was shown that elastic scattering cross sections are essentially 
independent of collision energy E at low energies, in electric fields sufficiently strong to 
polarize the molecules At the same time, state-changing cross sections scale as E^^"^ for 
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fermions and as for bosonic molecules. Therefore, at "sufficiently low" temperatures, 

elastic scattering is always larger for fermions, and evaporative cooling should be possible. 
Using the first Born approximation(BA), it was concluded that this is the case for the 
molecules OCS and CH3CI, at reasonable experimental temperatures. However, the BA 
may not be strictly applicable il8| for all fields and energies of interest. Indeed, the Fermi 
suppression of inelastic collisions may not be of great use for evaporatively cooling the OD 
radical Q- 

We are therefore motivated in this paper to re-visit the question of field- dependent scat- 
tering of fermionic molecules, from the perspective of close-coupling (CC) calculations. 
A complete theoretical description of molecule-molecule scattering is complicated by the 
complexity of the short-range interaction between molecules. This interaction is generally 
unknown to sufficient accuracy for cold collisions. Therefore in order to avoid the inclusion 
of unknown parameters of interaction, we seek and explore situations in which the infiuence 
of short-range physics is minimal. It appears that for weak- field seeking states the influ- 
ence of the short-range potential is suppressed, owing to avoided crossings in the long-range 
interaction For collisions of identical fermionic molecules, the infiuence of short-range 
physics may be even smaller, since only partial waves with / > 1 are present, and there is a 
centrifugal repulsion in all scattering channels. 

In this paper we seek a case where the elastic rate can be significantly larger than the 



me 



astic rate, by at least two orders of magnitude, and so provide feasible evaporative cooling 
21|. The most important finding here is that a large ratio of elastic and inelastic rates can 
be found higher temperatures where the Wigner threshold law is not valid. Theoretically, we 
suggest that this phenomenon is quite general, although large electric fields may be required. 
This effect exists even at temperatures around ImK which is already experimentally attained 
^|. It is in just such a gas that evaporative cooling would be a useful technique for achieving 
even higher phase-space density. 
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II. MODEL 



A. Polar ^S- type molecules 



The majority of diatomic molecules have electronic ground states 22|. The energy 
levels of these species can be described by the rotation J, total spin F (i.e., including 
nuclear spin), and vibration v quantum numbers. In this paper for simplicity we will neglect 
hyperfine splitting as the hyperfine interaction for molecules is smaller than for 11 or '^S 
molecules and we consider them only in the t; = vibrational ground state. So we will treat 
polar molecules as rigid rotors with a permanent dipole moment. The Stark splitting will be 
characterized by (J, Mj), where Mj is the projection of J on the direction of the external 
electric field. Thus the Hamiltonian for a polar molecule in a field is 

H ^ = Hrot + H field (1) 

The matrix elements for the Hamiltonian (Q) in this basis are 

< JM,|if^^|J'M,, >=i?eJ(J + l)5jj'-/i^(-l)"'-^([J][J'])^/M ^ ^ ^ I ^ ^ ^' (2) 

\^ ^ \^ -Mj Mj j 

In this expression is the rotational constant, ^ is the molecular dipole moment, £ is the 
strength of the electric field 

The different values of the molecular rotation J are strongly mixed in laboratory strength 
fields. Accordingly, in practice we transform the molecular state to a field-dressed basis for 
performing scattering calculations: 

\{J)Mj-£>=Y,a{J)\JMj>, (3) 
J 

where a (J) stands for eigenf unctions of the Hamiltonian determined numerically at each 
value of the field. The J quantum number is not a good quantum number in a field, but we 
will continue to refer to molecular states with J as a reminder of the zero-field value of J. 

Figure 1 shows the Stark energies computed using all the ingredients described above. In 
zero field the energy levels are determined by the rotational constant 3^.. We demonstrate the 
Stark splitting for linear OCS molecule which is quite typical for this type of molecule. For 
the other molecules we consider, the figure would look exactly the same, but with rescaled 
axes. We estimated that the hyperfine splitting is of order of ixK, which means that this 



effect may be important for ultracold energies. But here we will ignore the hyperfine effects 
as the region of energies we are considering is quite above ^K. The Stark shift is quadratic 
for fields below the "critical field" defined by £q = BeJ{J +l)/2fi. It is a rather approximate 
estimate because of the mixing between the neighboring rotation levels. For a molecule in 
its lowest weak- field -seeking state |10 > the critical electric field is typically on the order of 
10^ — 10^ V/cm for the species we consider here. For fields larger than this, the states with 
J = 0, 1, 2 are deeply mixed. As a consequence, a state like |10), which is weak- field seeking 
in low fields, can become high-field seeking at somewhat higher fields. In the following, as 
we are interested in J = 1 states, the critical field is given by Sq = Be/ fj,. 



B. Dipole-dipole interaction 

The intermolecular dipole- dipole interaction has the form: 



i?3 

= -^E(-irc'2.HK®H'- (4) 

where loi^2 = {^1,2, ^1,2) are the polar angles of molecules 1 and 2 with respect to the lab-fixed 
quantization axis, and R = {R, u) is the vector between the center of mass of the molecules 
in the laboratory fixed coordinate frame. Here C^g(co') is a reduced spherical harmonic jigj ]. 
We express the Hamiltonian in a basis of projection of total angular momentum, 

M = Mj,+Mj,+Mr, (5) 

Ml is the projection of the partial wave quantum number / on the laboratory axis. In this 
basis the wave function for two molecules is described as: 

^^ = 4 E {\1> ®\2> ®\lMi>}^ x^^'^'^R), (6) 

where {...j"^ is the angular momentum part of this wave function and \i > is the wave func- 
tion for each molecule. As we consider the target and the projectile as identical molecules, 
we must take into account the symmetry of the wave function (jU)) under exchange. 

Taking into account the Wigner-Eckart theorem, we can present the reduced angular 
matrix element as 

<121Mi\\Aa\\1'2'1'Mi, >= (7) 



M 



J2 



J[ 1 J I 



Mv 




1 2 

- Mj/ Ml - Ml, 

1 Ji J[ 

Mj, - Mj. -Mj, Mj, 

V L I 

Ml, Ml - Ml, -Ml 



X 



^ I' 2 I 







In practice, before each scattering calculation the Hamiltonian matrix has to be transformed 
from this basis into the field-dressed basis defined by Q. We solve the coupled- channel 



24| to calculate total state-to- 



equations using a logarithmic derivative propagator method 
state cross sections. Since the projection of total angular momentum on the field axis, Ai, 
is a conserved quantity, calculations can be performed for each value of M. separately. We 
find, generally, that the dominant contribution to cross sections arises from the minimum 
allowed value of Ai and that the general behavior of cross sections for other Ai is quite 
similar, and so restrict calculations accordingly. 

The scattering calculations quickly become computationally expensive as more rotational 
states and partial waves are included. Because we need to calculate cross sections at many 
electric field values, we choose a "compromise" basis set that includes rotational levels up 
to Jmax = 3 and partial waves up to Lmax = 3. This basis set then consists of 182 scattering 
channels. This level of approximation tends to get the general magnitude of elastic scattering 
cross sections fairly accurately, and to overestimate inelastic scattering cross sections. We 
therefore expect to draw conservative conclusions on the high ratio of elastic to inelastic 
scattering. 



III. RESULTS AND DISCUSSION 



We have chosen a variety of different linear molecules for this study, to span a ran ge o f 
rotational constants and dipole moments, and also to connect with the results of Ref . j20| . 



Properties of the molecules are summarized in Table I [25|, |26|. Of particular interest is 
the alkali dimer RbCs, which is a leading candidate in the experimental quest to observe 
cold collisions [2^. For these molecules, we are interested in the lowest energy weak-field- 
seeking state of the ground vibrational state, \JMj >= |10 > (see Fig^). In contrast to 
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29|, here it is necessary to include several 



the n molecules we have previously studied [28, 
rotational levels of the molecule. 

The general behavior of the thermally-averaged cross sections versus temperature is shown 
in Fig. 121 This example is for CICN molecules in an electric field of £^ = 20 kV/cm, well 
above the critical field for this molecule, meaning that the molecule is strongly polarized. 
The heavy solid and dashed lines represent (respectively) the elastic and state-changing 
cross sections for fermionic isotopomers of this molecule. It can be seen that the standard 
threshold behavior of colliding dipoles occurs at energies below several tens of microKelvin. 
Namely, the elastic cross section becomes a constant, and the inelastic cross section goes to 
zero as By contrast, the light solid and dashed curves show the same quantities, 



but for a bosonic isotopomer. Here the threshold laws work against the experimentalist, with 
the inelastic rate diverging as At higher temperatures, above the threshold 

regime, the elastic and inelastic cross sections are comparable for both bosons and fermions, 
consistent with what was found for OH and OD radicals in Ref . Q| . 

Remarkably, this situation can change quite dramatically if the electric field takes certain 
values. We illustrate this in Fig. 13 by plotting elastic and inelastic cross sections at T = 1 
/iK, versus electric field. In each case, the inelastic cross sections (dashed lines) fall by 
2-4 orders of magnitude at a specific value of the electric field. Here we show only results 
for fermionic species. For bosonic species the corresponding drop is smaller because of the 
threshold behavior at this energy. The field values at which this suppression occurs, denoted 
£^supp, are tabulated in Table I. In all cases the suppression occurs at around 3 times the 
critical field So- To explore this phenomenon further, we focus on a particular example in 
the following, namely, a fermionic isotopomer of CICN. 

The unexpected suppression of inelastic collision rates has its origin in the interplay be- 
tween the strong dipole-dipole interaction and the rotational energy levels of the molecules. 
The first clue to the mechanism of suppression comes from considering the scattering thresh- 
olds, i.e., the total energy of both molecules when they are far apart. (Fig. |^. Here 
each threshold is labelled according to the quantum numbers | JiMi, J2M2) of the colliding 
molecules, and zero energy represents the lowest energy threshold corresponding to |00, 00). 
The solid line denotes the |10, 10) incident channel of interest to this paper. It can be seen 
here that this threshold crosses three other thresholds, for the channels |00, 22), |00, 21) and 
|00,20), at fields S = 11.6,12.8,13.65 kV/cm respectively, very near to £^supp As we will 



7 



see below £^supp can be defined more quantitatively as the field at which thresholds for the 
channels |10, 10) and |00,20) are degenerate. 

This is a strange situation, in which the state-changing collision cross sections diminish 
sharply just as new states become energetically available. The second piece to this puzzle is 
found by examining the approximate adiabatic potential energy curves for two situations, as 
shown in Fig. |S| To simplify these curves, we have included only the partial wave L = 1 in 
their construction, although L = 3 is also used in the multichannel scattering calculations. 
In (a) are shown the curves for S = 12.3 kV/cm, just below the threshold crossing. Here the 
incident channel |10, 10) is below the nearby thresholds, and correlates adiabatically to the 
solid curve. This curve draws molecules into the small- i? region where they interact strongly 
and can readily change their internal state. 

By contrast. Fig. Efb) shows the adiabatic curves at a field S = 13.7 kV/cm, just above 
the threshold crossing. Now the incident channel correlates adiabatically to a repulsive curve, 
so that the molecules do not approach one another nearly as closely as in the previous case. 
This shielding, in turn, reduces the likelihood of inelastic collisions. The dominance of the 
repulsive curve is a peculiarity of the strong, anisotropic dipole-dipole interaction, and is 
similar to the physics that generates "field-linked" molecular dimer states, which also keep 
the molecules far from each other 



28[. While this simple picture of the suppression is 
probably a good first approximation, there is clearly more going on. This can be seen in the 
tabulated values of £^supp, which are not always equal to 3£o, even though this is the field 
value where the thresholds cross for any molecule. 

We also remark that we have noted a similar, but somewhat less dramatic, suppression of 
inelastic cross sections in bosonic analogues of the molecules considered (Fig. Presumably 
fermions have an advantage since they have nonzero partial wave angular momentum in all 
channels, which aids in the shielding effect discussed above. 

The strong suppression in the low-energy limit naturally has consequences at higher 
collision energies. Figure U\ shows the elastic (solid) and inelastic (dashed) cross sections 
versus temperature, at a field S = 14 kV/cm, where the inelastic rates have just become 
suppressed. Strikingly, the ratio of elastic to inelastic cross sections is close to two orders of 
magnitude even at temperatures as high as several mK, which is easily attainable in Stark 
deceleration or buffer-gas cooling experiments. For this reason, it is conceivable that cold, 
dense samples might be amenable to evaporative cooling that will reduce them to ultracold 
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temperatures. 

Figure [7| also illustrates the effect of increasing the number of channels in the scattering 
calculations. For both elastic (solid lines) and inelastic (dashed lines) cross sections, four 
curves are shown, corresponding to various maximal numbers of rotational states (Jmax) 
and partial waves (Lmax)- In all cases, changing the size of the basis set has little influence 
on the overall elastic scattering cross section, although the features understandably shift in 
field. By contrast. The inelastic cross sections are quite sensitive to Jmax, dropping more 
than an order of magnitude as Jmax increases from 3 to 4. In addition, the number of partial 
waves plays a role in the actual cross sections. In any event, the conclusions still hold, and 
the ratio of elastic to inelastic cross sections should be quite high. 



IV. CONCLUSION 

In this paper we considered the coUisional dynamics of polar molecules in a dc-electric 
field taking CICN, HCN, OCS and RbCs as prototypes. The lowest weak-field-seeking 
state was studied from the point of view of evaporative cooling. As a rule the strong and 
anisotropic dipole-dipole interaction should provide quite large inelastic rates. We have 
found, however, that polar state molecules possess a substantial suppression of inelastic 
collisions at high values of the applied electric field. The sudden drop in inelastic cross section 
coincides with the degeneracy of certain molecular rotational levels. Adiabatic pictures 
reveal that the interaction changes from mostly attractive to mostly repulsive upon crossing 
the field where this coincidence occurs. The strong suppression of inelastic scattering from 
the \JMj) = 1 10) state of a molecule seems generally to occur at a field nearly equal to 
3i?e/ /i, and may enable evaporative cooling in an electrostatic trap. 

As a final remark, we note that in many experiments alkali dimer molecules are formed 
in high-lying vibrationally excited states, as a result of photoassociation or magneto- 
association. In these quite different states, the field scale can be significantly larger. To make 
an estimate, consider the RbK molecule, which is being pursued by the UConn group js^. 
In a hypothetical weakly-bound state whose outer turning point is 40 a.u., this molecule's 
rotational constant is approximately = 1.4 cm^^. The dipole moment of this molecule 
has been estimated to be approximately 1.4 x 10^^ eao at this intermolecular separation jsi], 
yielding an critical field of ~ 2 MV/cm. Thus it would seem unlikely that the effects de- 
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scribed here are observable for such weakly-bound molecules. For even more weakly-bound 
molecules, the critical field quickly becomes larger, owing to the exponential falloff of the 
dipole moment. 

This work was supported by the NSF and by a grant from the W. M. Keck Foundation. 
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TABLE I: Molecular parameters for the species considered. The rotational constants B^, and dipole 
moments fi for the triatomics come from Ref. [28], while those for RbCs come from Rcf. [29]. Here 
So = Be/fi is the "critical" field, while £"supp is the calculated field value at which the inelastic 
processes become suppressed. 

Molecule Be (cm-^) /x (D) SSq (kV/cm) S^npp 

RbCs 0.017 1.26 2.35 2.55 

CICN 0.199 2.833 12.57 13.65 

OCS 0.203 0.715 50.71 55.3 

HON 1.478 2.985 88.51 96.4 
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E(K) 

FIG. 2: Elastic(solid lines) and inelastic(dashed lines) cross sections for CICN molecule at an elec- 
trostatic field E = 20000y/ cm. Thick and thin curves are for Fermi and Bose particles respectively. 
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FIG. 3: Elastic(solid lines) and inelastic(dashed lines) cross sections versus an electrostatic field 
at collision energy for a)ClCN, b)OCS, c)HCN and d)RbCs molecules. For these calculations, 

Lmax — 3, Jjnax — 3 
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FIG. 4: Threshold energies for CICN molecules referred to the threshold of the |00, 00 > channel. 
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FIG. 5: Adiabatic curves for CICN molecule at a)12.3 kV/cm and b)13.7 kV/cm. In both cases, 
the dark curve correlates adiabatically to the |10)|10) incident channel at large K. 
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FIG. 6: Elastic(solid lines) and inelastic (dashed lines) cross sections versus electrostatic field at 
collision energy IfiK for CICN molecule. Thick and thin curves are for Fermi and Bose particles 
respectively. Mtot = 
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FIG. 7: Elastic(solid lines) and inelastic (dashed lines) cross sections versus the collisional energy 
for CICN molecule. Jmax and L^ax are maximal taken into account values of rotation an partial 
quantum numbers. The electrostatic field is \AkV/cm 
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